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Abstract 

Using Dyson-Schwinger equations within an approach developed by 
Broadhurst and Kreimer and the renormalization group, we show how 
high loop order of the renormalization group coefficients can be efficiently 
computed in a supersymmetric model. 

1 Introduction 

Perturbative Quantum Field Theory is known for its tremendous successes, 
with its ability to obtain highly precise values in Quantum Electrodynamics, or 
the tests of the standard model involving radiative corrections to weak interac- 
tions. However, actual calculations become rapidly cumbersome and display a 
conjunction of analytical and combinatorial difficulties combined with the very 
fast growth of the number of relevant terms, which get compounded by the need 
to address the subtraction of the various sub divergencies. 

In the quest of organizing principles for taming the combinatorial problem, 
a major progress has been the recognition of a Hopf algebra structure in the 
renormalization of Quantum Field Theories pQ and subsequent developments 
2J- 4J and the link relating the cohomology of the introduced Hopf algebras with 
Dyson- Schwinger equations [5]- [6] (See [7] for a complete list of references). An 
application to the summation of a category of diagrams in a simple Yukawa 
field theory or a 3 theory in 6 dimension has been obtained by Broadhurst 
and Kreimer in [6]. In that work a very powerful method to compute higher 
order corrections to the anomalous dimension has been developed. Basically, it 
is based on a differential equation derived from the Dyson-Schwinger equation 
for the self-energy. 
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The purpose of this letter is to extend such results to the case of a super- 
symmetric model. We have considered a Wess-Zumino like model in which 
the non-renormalization of the vertex ensures that the order zero in a large N 
expansion of the renormalization group functions can be obtained. We also ex- 
pect that the known non-renormalization properties related to supersymmetry 
should allow to go beyond this initial category of diagrams in future works. 
Along the way, we propose an alternative derivation of the propagator-coupling 
duality presented in [6] based on intuitive renormalization group arguments. 



2 The model 

We consider massless chiral superfields ^ and (i = 1,2,... ,N) and their 
(antichiral) complex conjugates and which satisfy the following con- 
straints: 

D a * + = 0, D a $f = (1) 

The derivatives are defined as: 

in terms of the chiral coordinates y defined as 

yf = x ^ + 19^9 (3) 

Each chiral superfield represents a complex scalar (A,Bi), a Weyl fermion 
(Xi £i) an d a complex auxiliary field (F, Gi) as it can be seen from their expan- 
sion in the 9 variables, 

* = A(y) + V29 X (y) + eeF(y) 

*i = Bi{y) + V26Uv) + QOG^y) (4) 

Dynamics is governed by the Lagrangian density 

„ N . N . 

L= d 4 9W + + J2 d 4 9$i<S>t + -^=Y J d 2 9^ 2 (5) 
i=l ^ i=1 J 

The cubic superfield interaction corresponds to a Yukawa coupling between 
scalar and fermions and, after the elimination of the auxiliary fields, to a quartic 
coupling of the scalars. The Lagrangian reads in terms of component fields: 

N N N 

i=l i=l i=l 

N 

+J= Y {A 2 Gi + lABiF - X 2 Bi - 2^ iX A + h.c.) (6) 

* i=l 

In our conventions (the same as in [8]), we have diag g mn = (—1, 1, 1, 1) and 

□ = -d 2 t + V 2 (7) 
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Figure 1: Schwingcr-Dyson equations 

3 Dyson— Schwinger equations 

We only consider corrections to the propagators of the fields associated to the 
superfield as corrections to the propagation of the superfield <&j get a factor 
1/N and do not contribute in the large N limit. Moreover super symmetry 
ensures that there are no divergent vertex corrections so that what we are 
calculating represents effectively the order zero of the large N development. 

The set of diagrams that we consider is the same as the one studied in [BJ, 
except that now the diagrams have to be considered as superfield diagrams and 
therefore, each one will correspond to a sum of ordinary Feynman diagrams. 
The number of such diagrams grows rapidly with the number of loops and for 
example a given 12 loop superdiagram corresponds to 18144 ordinary diagrams. 

The diagrams we are looking at correspond to the expansion of the Dyson- 
Schwinger equations graphically represented in Figure 1. 

With our (Minkowski) conventions, free propagators read 

n-» = (0\T(A(x)A*(x'))\0)= l n~ 1 (x-x') 
D^(p) = (0\T(F(x)F*(x'))\0) = i5(x-x') 

(V£)jp) = (0\T{ Xa (x)x $ (x'))\0)=a^d m D~ 1 (x-x') (8) 

It will be convenient to write Schwinger-Dyson equations in Euclidean space. 
Then calling Up, H x , Ha the full propagators for fields F, x and A in (Euclidean) 
momentum space, the corresponding self-energies are introduced as 

u- A Hp) = p 2 (l-S A (p 2 )) 

yt f \ p ) = -(i-s F ( P 2 )) 

n-Hp) = Pm a m {i-j: x (p 2 )) (9) 

We consider a simple one-loop correction to the self-energy, where one of the 
propagators has a momentum-dependent renormalization T,(p 2 ). The simplest 
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case is the one for the self-energy of the auxiliary field F, since there are only 
contributions from a single diagram with scalar propagators: 

W) = -^/^V(i-s^))(p-^ (10) 

For the fermion Xi we have the two contributions shown in figure 1, 

n 2 / f n rr m 



Pma Sx(p) = [J dq w^wm^w + 

V(i-X4(g 2 ))(p-g)V 1 j 

In the case of the scalar ^4, the three different terms shown in figure 1 evaluate 
as: 



47r 4 V7 V(i-£a(<? 2 )) '7 "(i-sh^Kp-?) 2 

, f a -Tr(q n a n (q m a m -p m a™)) \ 

+ ] dq q 2(i-z x ( q *))(p- q y ) [U) 

A consistent ansatz to solve these equations is to take 

Ha = S x = Tip = X (13) 

Indeed, with this ansatz, the two integrals in (jlip combine to give Pmcr m times 
the right hand side of (|10p . Similarly, the equation for (|12p becomes p 2 
times ([TO]) , Note that Tr(cr n cr m ) = —2g mn . 

The integration in (jlOp can be done along the same lines as in [6]. The 
angular integration uses the fact that the angular average of l/(p — q) 2 is 
l/max(p 2 , g 2 ). With the variables x = p 2 and y = q 2 , the radiative correc- 
tion becomes 



1 Z 100 1 



(14) 



/ a x(l-X(y)) 7 2 y(l-E(y)) 

Taking the x derivative to get rid of the infinite constant, we arrive first at 

d a 2 f x 1 

x 2 — E(x) = / dy - - (15) 



w 4vr2 7 "(l-E(y)) 
A further derivation allows to obtain a simple differential equation for X, 

2 

(1 - S(x))Z>(Z) + l)E(z) = ^ = a ( 16 ) 
with D = xd/dx, which is equivalent to D = d/d log(p 2 / fi 2 ) . 
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4 Anomalous dimension 



The equation (|16p is very similar to the one obtained for the Yukawa case in |6j . 
The only difference is that the operator D{D + 1) replaces in the present case 
the D{D + 2) one in [6j. The two equations can then be identified by rescaling 
the variable L = log(p 2 //i 2 ) and a by appropriate factors. The results derived 
in [6] can then be reobtained taking into account the correct factors. 

The most direct way of obtaining a solution of eq. (|16p is to solve for 
order by order in a, using the renormalization constraint S(/x 2 ) = 0. At order 
n in o, Ti(x) is a polynomial in L of order n without constant terms. The 
renormalization group 7 function is then the first derivative in L for L = 0. 



7(«) 



dlog(l-£tp 2 ,a)) 



dlog(p 2 ) 



(17) 



Here 7 is the anomalous dimension for all the fields in the scalar multiplet ^. 

As was shown already in [6\, this procedure can be further simplified if 
one directly computes 7 through the use of a non-linear equation that can be 
deduced from eq. (|16p . Indeed, consider the renormalization group relation [9] 



d 



d{^,a(s 2 )^d(^. 



where 



d{p 2 , a)' 1 = 1 -S(p 2 ,a) 
Differentiation of eq. (|18p with respect to p 2 one obtains 



dlog(l 



or 



d log p 2 
dlog(l-E) 



7(a(s 2 )) 
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(18) 



(19) 



(20) 



(21} 



dlogp 2 1 \ (1 

This last identity was obtained as a consequence of the non-renormalization of 
the vertex which implies a connection between the (3 and 7 functions of the 
renormalization group. It should be noted that eq. (|21|) . which was derived 
here using the renormalization group, can also be derived using Hopf algebra 
arguments as explained in [6]. Indeed, in this last work the dependence of the 
self-energy on the momentum and the anomalous dimension (the propagator- 
coupling duality) based on the Hopf structure of Feynman diagrams is used to 
prove eq. (|2ip . Here, our derivation relies instead on the renormalization group 
relation eq. (fl8l) . 

Let us now rewrite eq. f)21 1) in the form 



D(l-X) = (l-£) 7 (a/(l-£) 2 ) 



(22) 



and apply the operator {D + 1) to both sides. After some work one ends with 
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-a + a 



da 



(23) 
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From this nonlinear differential equation one can derive the perturbative ex- 
pansion of 7, 
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n>0 
n 



Gn+1 = GfcGn+l-fc (24) 

k=l 

A simple Mathematica program allows to evaluate G n for large values of n 
according to this recurrence relation. As an example, Gioo can be evaluated 
quasi immediately giving 

G n = 2421240422948166903917318482712165846373737073005111621110 
3204862084235299310224396534696102431673533019446585512533 
0800645887484858421550031445472180813110312557469722336097 
3401051786164 (25) 

We can then compute 7 as a function of a using eq. (|24p . However being 
G n ~ (2re — 1)!! for large re, the series is asymptotic and therefore we use a 
Pade-Borel resummation method. We start from the approximation 

7 — a J dxP(ax)\/ r x exp(-x) (26) 

where P(x) is defined by the following series 

^>=£(-D"f^) (27) 

Numerical values for 7 can then be obtained using a Pade [iV\M] approximant 
for P(x). If one avoids the values of [iV\M] for which this approximant has 
poles on the positive real axis, one can obtain stable evaluations for growing 
values of a. The successive approximations are fully coherent and higher degree 
approximants allow to reach higher values of the coupling constant g, a = 
g 2 /Air 2 . In a few seconds one can compute the anomalous dimension 7 for 
values g 2 ~ 700 with a relative accuracy of the order of 10 -10 . As an example, 
one finds ^(g 2 = 789.568) = —6.30706 with a [99\100] approximant. Rational 
solutions are not possible for such high orders and the precision afforded by 
floating point numbers is not sufficient to obtain solutions, but the use of higher 
precision allowed to reach easily this stage. 

In figure 2, we show 7 as a function of g 2 . Remarkably, we see from the 
numerical results that the anomalous dimension 7 of the scalar field can reach 
values where it exceeds the canonical one, so that, unexpectedly the propagator 
grows with the impulsion. One confirms this analytically from the asymptotic 
behavior of 7 for large g that can be easily inferred from eq. (|23p . 

7 « gy^ (28) 

We see that our calculation gives a complete knowledge of the renormal- 
ization group functions since, due to the non-renormalization of the vertex one 
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-y 




has (5 = —2"f. In particular, the asymptotic behavior of 7 in eq. ()28[) allows to 
prove that the coupling constant goes to infinity for a finite value of the scale. 

We end this section by noting that the renormalization group played a cen- 
tral role in obtaining the basic equations (I2ip - (|23p from the Dyson-Schwinger 
equations. General information on the renormalization group have been pre- 
sented in [4] , deduced from the Hopf algebra structure of perturbative quantum 
field theory. It would be of interest to combine the two approaches to fully real- 
ize the potential of Dyson-Schwinger methods to obtain renormalization group 
functions. 



5 Discussion 

In this letter we have generalized the calculations of [6] to a supersymmetric 
setting. We have chosen to study a Wess-Zumino like model in which the 
vertex is not renormalized and this allowed us to calculate the order zero of the 
renormalization group functions in a large N expansion in a closed form. 

Writing the Dyson-Schwinger equations for propagators and proposing an 
appropriate ansatz for the corresponding self-energies we have found a non- 
linear differential equation for the anomalous dimension (eq. (|23p ) from which 
one can derive a perturbative expansion for 7. Coefficients in such expansion can 
be very easily evaluated up to very large orders and a Pade-Borel resummation 
allowed us to compute 7(5) in a large range of coupling constants with great 
accuracy. 

In our calculations, the power of the Hopf algebra approach has not be used 
explicitly. It is tempting to speculate, however, that the special features of 
perturbation theory in supersymmetric models should reflect in some particular 
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structure within the Hopf algebra approach. We intend to discuss this issue in 
a future work. 
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